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Abstract 

It is well known that when the distance t traveled by neutrinos in matter is short, 
matter effects in oscillations of neutrino flavour are small and decrease with decreasing 
t more rapidly than the oscillations effects themselves. We discuss the reason for 
this and demonstrate that under certain circumstances this statement is no longer 
correct. In particular, we show that if neutrinos propagate significant distances in 
vacuum before entering matter (or after exiting it), matter effects in short-i neutrino 
oscillations can be significantly enhanced. Implications for oscillations of solar and 
atmospheric neutrinos with nearly horizontal trajectories inside the earth and for 
accelerator experiments are considered. We also comment on neutrino oscillations in 
matter due to flavour changing neutral currents. 
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1 Introduction 



It is well known that when the distance t traveled by neutrinos in matter is short, the 
probabilities of oscillations of neutrino flavour in matter reduce to those in vacuum, i.e. 
with decreasing t matter effects on neutrino oscillations die out more rapidly than the 
oscillations effects themselves. In the present Letter we show that under certain conditions 
this is no longer true and matter effects in short-t experiments can be quite significant or 
even dominate the oscillation probability. In particular, we show that if neutrinos propagate 
significant distances in vacuum before entering matter (or after exiting it) matter effects 
in short-t neutrino oscillations can be strongly enhanced. We also discuss oscillations of 
neutrino mass eigenstates in matter, which are relevant for oscillations of solar and supernova 
neutrinos inside the earth. We show that these oscillations can have sizeable probabilities 
even when neutrino pathlengths in matter are relatively short. 



2 Two-flavour neutrino evolution in matter 

Consider the evolution equation for two-flavour neutrino oscillations in matter in the weak 
eigenstate basis [|I|, [|: 

4UH-r ^>)(:H<>(:;) « 

Here v a ^{t) are the probability amplitudes of finding neutrinos of the corresponding flavor 
a, b at a time t (in particular, one of these two species can be a sterile neutrino v s ). The 
parameters A and B are 

A(t) = cos 20 o 6 - V (t) , B = sin 29 6, (2) 

where 

*=^- nt) = %N(t). (3) 

Here Gf is the Fermi constant, E is neutrino energy, Am 2 = m 2 , — m\, where are the 
neutrino mass eigenvalues, and 9$ is the mixing angle in vacuum. The effective density N(t) 
depends on the type of the neutrinos taking part in the oscillations. For transitions between 
antineutrinos one should substitute —N for iV in eq. (0). In a matter of constant density 
the probability of v a <-> oscillations takes a very simple form 

P(v a -> Vb\ t) = sin 2 29 sin 2 ut , (4) 

where 9 and 2uj are the mixing angle in matter and the energy splitting between the matter 
eigenstates respectively: 
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sin29 = — = sin2# -, lu = VA^+W . (5) 
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It is now easy to see that in the short baseline limit (cut <C 1) the oscillation probability 

/\ 777 2 

P[v a -> u b ; t) ~ sin 2 26 (cut) 2 = sin 2 2# {6 -tf = sin 2 2# (^^) 2 , (6) 

which is just the short-baseline limit of the oscillation probability in vacuum [[|. 
This raises a number of questions: 

• Why do the matter effects on oscillations die out with decreasing t more rapidly than 
the oscillation effects themselves? 

• Is this true also in a matter of non-constant density? 

• Is this true in the case of oscillations between more than two neutrino species? 

• Are there any conceivable situations when the matter effects on oscillations of neutrinos 
traveling short distances in matter can be as large as the oscillation effects themselves? 

The answer to the first of these questions will help us to find the answers to the rest of 
them. 



3 Short-pathlength neutrino oscillations in matter 

Consider the two-flavour evolution equation ([]]) in the short baseline (small oscillation phase) 
limit. In this limit the oscillation effects are known to be small, so one can solve the evolution 
equation in perturbation theory. Assume that the initial neutrino state is v a . Then in the 
leading order in perturbation theory the amplitude Vbif) — —iBt, which immediately yields 
Thus, the fact that the matter effects on neutrino flavour oscillations disappear in the 
short baseline limit is the consequence of the two facts: (1) the oscillations effects are small 
and the perturbation theory applies; (2) since the initial state is a flavour eigenstate, in 
the leading order the transition amplitude is only determined by the off-diagonal terms in 
the effective Hamiltonian in (]T|), whereas the matter effects enter only through the diagonal 
terms. 

From the above the answers to the other questions that we asked immediately follow: 
The functional t dependence of matter density is irrelevant, i.e. our conclusion holds ir- 
respective of whether the matter density is constant or not; it is also true in the case of 
oscillations between more than two species because, in the weak eigenstate basis, matter 
density enters only into the diagonal terms of the effective Hamiltonian for any number of 
flavours; the above conclusions need not be correct if the initial neutrino state is not a weak 
eigenstate, i.e. not a neutrino of a definite flavour. 

Let us now consider the situation when the initial neutrino state entering the matter is 
a coherent superposition of the flavour eigenstates u a and v^. \ui) = \v(Q)) = a\\u a ) H-g^z^), 
I a-i 1 2 + 1 02] 2 = 1. In what follows, unless otherwise is specified, we shall always understand 
by neutrino pathlength the distance traveled by neutrinos in matter. It is straightforward 
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to find the transition probabilities to the second order in perturbation theory: 



Pivi -> v a \ t) ~ |ai| 2 -(|ai| 2 -|a 2 | 2 )( J Bt) 2 -2Im(a;ai)( J Bt)-4Re(a;a 1 ) J B / dh [ * dt 2 A{t 2 ) , 

Jo Jo 

(7) 

P(i/i -> f 6 ; t) ~ |a 2 | 2 + (|ai| 2 -|a 2 | 2 )( J Bt) 2 + 2Im(a;ai)( J Bt)+4Re(a;a 1 ) J B / tftj / * dt 2 A(t 2 ) , 

Jo Jo 

(8) 

For our further discussion it is useful to write down the simplified versions of these expres- 
sions relevant for neutrino oscillations in a matter of constant density: 

P(u { -> v a \ t) ~ | ai | 2 - (|ax| 2 - \a 2 \ 2 )(Bt) 2 - 21m(a* 2 a 1 ){Bt) - 2Re(a* 2ai ) (AB t 2 ) , (9) 

P(ui -> i/ 6 ; t) ~ |a 2 | 2 + (|ax| 2 - \a 2 \ 2 )(Btf + 2\m(ala l )(Bt) + 2Re{a* 2 a l )(AB t 2 ) . (10) 

Matter effects enter through the parameters A in the last terms on the r.h.s. of eqs. (|7|)- 
(|10|). Notice that the leading order matter-induced contribution is ~ V5-t 2 , whereas in the 
case when the initial state is a pure flavour eigenstate the leading order matter effect is 
~ (V5-t 2 ) 2 (assuming V > 5). Thus, matter effects in short-pathlength neutrino oscillations 
can be strongly enhanced when the initial state is not a pure flavour eigenstate []. 

As can be seen from (|7|)- ([T0|) , the oscillation effects themselves get enhanced in this case 
provided that Im^ai) 7^ 0, the leading contribution to the transition probability being 
now ~ S-t rather than ~ {5-t) 2 . 



4 Two-media neutrino oscillations 

How can one create an initial neutrino state which is not a pure flavour eigenstate? One 
possibility is to let neutrinos oscillate in a different medium before entering the medium 
of interest. For example, if neutrinos propagate in vacuum before they enter the matter, 
the initial state arriving at the vacuum-matter border is no longer a flavour eigenstate but 
rather a coherent superposition of the flavour eigenstates. 

Assume that neutrinos are initially produced in the flavour eigenstate u a , propagate a 
distance t\ in vacuum and then the resulting state propagates a distance t in a matter of 
constant density. The oscillation phase acquired in vacuum is <p\ = 8-t±. Let us denote 
sin 0i = s±, cos 0i = c\. Then the state |z/j) entering the matter is characterized by a\ = 
c\ + i cos 2#o Si , a 2 = —isin26 si. In the limit of short-pathlength oscillations in matter 
one finds from ([]) and (|TUD Q 

P(ui -> v h \ t) ~ sin 2 26 [s\ + (c 2 - s\)(5-t) 2 + 2s l c l {5-t) + 2s\ cos 29 (V5-t 2 )} (11) 

1 It is also important that the term VS-t 2 is sensitive to the sign of Am 2 , whereas (VS-t 2 ) 2 is not. 
2 These probabilities can also be directly obtained (as a small fa limit) from the general expression for 
the evolution matrix for neutrino oscillations in two layers of different constant densities derived in B . 
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and P{vi — > z/ a ; t) = 1 — P(z/j — > t). Notice that the l.h.s. of eq. flTTD can also be 
understood as P(y a — > u b ; t\ +t). 

Eq. ([U]) has a simple physical interpretation. In the limit t — > one has P[i>i — > 14, ; t) = 
sin 2 2#o s 2 which corresponds to propagation only in vacuum; the second and the third terms 
in the curly brackets are due to the increase of the oscillation phase during the time interval 
t, neglecting the matter effects. These terms come from the expansion of sin 2 (0i + 5-t) in 
small 5-t. The fourth term is the leading order matter contribution. Notice that we do not 
assume the smallness of the phase 0i, and in fact a sizeable enhancement of matter effects 
is only possible when it is not small. 

Several comments are in order. 

(i) The leading-order matter contribution vanishes in the case of maximal mixing in 
vacuum, 6q = 7r/4. This may be useful for studying deviations of lepton mixing from the 
maximal one. 

(ii) The first three terms in the curly brackets in eq. are even in 5 while the last, 
matter- induced, term is odd in it and so is sensitive to the sign of Am 2 . [] Thus, enhanced 
matter effects can facilitate studying the type of the neutrino mass hierarchy. 

(iii) Although the absolute magnitude of matter effects, of course, strongly depends on 
8q, their relative size is only mildly 8q dependent provided that the vacuum mixing angle is 
not too close to 7r/4. This follows from the fact that sin 2 29o is the common factor in (|TT|). 

(iv) If the phase 0i = 5-t\ is not very close to 7r/2, the transition probability itself (i.e. 
neglecting matter effects) is also considerably enhanced as it now contains an ~ 5-t term in 
addition to the ~ (5-t) 2 ones. The relative size of the matter effects in this case is 

cos 29 tan fa (Vt) , (12) 

while for 0i « tt/2 it is 

cos26 (V/5). (13) 

Thus for generic values Si,Ci ~ 1 the contribution of the matter effects to the transition 
probability, though strongly enhanced compared to the case of pure flavour eigenstate en- 
tering the matter, is relatively small. It may, however, still be noticeable as the short 
pathlength approximation applies even for Vt as large as ~ 1/3 (the corrections are of the 
order (Vt) 2 ). When fa ~ 7r/2, matter effects give an important contribution to the transi- 
tion probability over the time period t; they dominate in the limit V 5. This fact can be 
used for experimental searches of matter effects in short-pathlength neutrino experiments. 

(v) The estimates ( |T2"D and ( |T3"D apply to the situations when the probability of finding 
Vf, at the vacuum-matter border is either experimentally known (two-detector experiments), 
or can be reliably estimated theoretically. If this is not the case, one has to compare the 
matter-induced contributions to P[yi h>b) with the probability of finding Uf, in the final 

3 More precisely, it is sensitive to the sign of cos 26*o Am 2 as only the sign of this quantity has a physical 
meaning. We adopt the convention cos 29q > and allow for both positive and negative signs of Am 2 . 
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state itself rather than with the increase of this probability due to neutrino propagation in 
matter. The relative contribution of matter effects is then 

cos2# (V<5-t 2 ). (14) 

For Vt ~ 5-t ~ 1/3 it can be about 10%. 

(vi) Eq. ([LTD is also valid when neutrinos first propagate a distance t in matter and 
then a distance t\ in vacuum. This is a consequence of the fact that two-flavour neutrino 
oscillations in matter are invariant under time reversal even if the matter density profile is 
not T invariant. Indeed, let the evolution matrix for eq. ([[]) be U(t2,ti), so that (^(£2)) = 
U(t2,ti)\v(ti)) . The time-reversed evolution matrix is Uifi^t^) = Ufojti)^ 1 = U(t 2 ,tiy. 
Since for any unitary 2x2 matrix | U21 \ = \ Uvi \ , the probability of neutrino flavour oscillations 
is T invariant. Notice that this is not in general true in the case of oscillations between n > 2 
neutrino species. 

Consider now a few numerical examples. 

Atmospheric neutrinos coming to a detector from below the horizon propagate first in 
the air (which for the purposes of neutrino oscillations can be considered as vacuum) and 
then in the matter of the earth. The distances that neutrinos travel in the atmosphere t\ 
and in the earth t are given by 

tt = -R\ cos 0| + y/(R + h) 2 - R 2 sin 2 6 , t = -2RcosB, (15) 

where R = 6371 km is the radius of the earth, is the zenith angle of the neutrino 
trajectory, and h ~ 15 km is the average hight at which the neutrinos are produced in the 
atmosphere. We are interested in the regime t\^> t which corresponds to nearly horizontal 
neutrino trajectories (zenith angles only slightly exceeding 7r/2). In this case one can have 
sizeable phases <p\ whereas neutrino oscillations in the earth are in the short pathlength 
regime. This corresponds to |cos0| < 0.01. Atmospheric neutrinos with nearly horizontal 
trajectories pass through the earth's crust where the density is nearly constant and equal to 
about 2.8 g/cm 3 and the electron number fraction Y e ~ 0.49 (the same is also true for short 
pathlength accelerator neutrino experiments). For oscillations between active neutrinos 
«-> P e or v e «-> v T this gives V ~ 5.17 x 10 -14 eV; for active-sterile neutrino oscillations V 
is a factor of two smaller. 

Consider, e.g., atmospheric «-» v e (or v e «-> u T ) oscillations for the neutrino trajectory 
with the zenith angle ~ 1.581 (cos© = —0.01). From eq. QT51 ) one finds t\ ~ 378 km, 
t ~ 127 km, which gives Vt ~ 3.3 x 10 -2 . For Am 2 = 3.2 x 10~ 3 eV 2 , which is the current 
best-fit value of the Super-Kamiokande atmospheric neutrino data |5[], and E = 1 GeV, 
one has 5 = 8 x 10~ 13 eV, <f>\ ~ 1.535, 5 ■ t ~ 0.517. The relative matter contribution to 
the total transition probability is 4.3% and that to the increase of the probability due to 
neutrino propagation in the earth is 18%. This is a large effect, taking into account that the 
distance neutrinos travel inside the earth is only 127 km. If neutrinos did not oscillate in the 
atmosphere before entering the earth, this contribution would have been only about 0.5%. 
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It should be noted, however, that the short-pathlength regime is valid only in a narrow 
range of zenith angles. In addition, the effective mixing angle for v e «-> v x oscillations is 
known to be small 0. 

Oscillations of into sterile neutrinos are disfavoured as the dominant channel of the 
atmospheric neutrino oscillations 0, but allowed as a subdominant channel with a weight 
that can be as large as about 50% |J. For this channel, relative matter effects are a factor 
of two smaller than they are for the <-> v e or v e <-> v T channels, but their absolute value 
can be significantly larger because the corresponding mixing angle can be quite large. 

Consider now the situation when neutrinos first propagate in matter and then in vacuum. 
This could, e.g., be realized in accelerator neutrino experiments with a detector placed 
on an earth's satellite, which is certainly a rather remote possibility. Let us assume for 
definiteness that the distance that neutrinos propagate in the earth t = 730 km (the baseline 
of CERN - Gran Sasso and Fermilab - Soudan mine experiments). This corresponds to 
cosG = —5.7344 x 10~ 2 , Vt ~ 0.191 (for oscillations between active neutrinos). Assuming 
that the height of the satellite's orbit is 750 km, the distance that neutrinos propagate after 
exiting the earth t\ ~ 2837 km. Let us again take Am 2 = 3.2 x 10~ 3 eV 2 . For E = 8 GeV, 
which is a typical energy of the accelerator neutrino experiments, one has 5 = 10~ 13 eV, 
<pi ~ 1.438, 5 ■ t ~ 0.37. The relative matter contribution to the total transition probability 
is 13%, and that to the transition probability acquired due to neutrino propagation in the 
earth is a factor 1.3. This means that matter effects dominate the transition inside the earth 
in this case. If the detector is placed on the surface of the earth, i.e. t\ = 0, matter effects 
constitute only about 3.5% of the oscillation probability for neutrinos that have traversed 
the earth, i.e. are almost a factor of 37 smaller. Thus, as paradoxical as it looks, the earth's 
matter effects on the oscillation probability become stronger when neutrinos are detected 
farther from the surface of the earth. 

Can one achieve a significant enhancement of the matter effects in the accelerator exper- 
iments by having very long decay tunnels? Unfortunately, this does not seem to be possible 
since a sizeable enhancement is only achieved when the distance traveled by neutrinos in 
vacuum is much larger than their pathlength in matter, and the latter should be at least a 
few hundred km. 

5 Day-night effect in solar neutrino experiments 

Another example of an initial state which is not a flavour eigenstate is given by the earth 
matter effects on solar neutrinos coming to a detector during night in the case of the MSW 
||T], Q solutions of the solar neutrino problem. In this case the neutrino state arriving at the 
earth is an incoherent superposition of the mass eigenstates V\ and u 2 (see, e.g., [jj for a 
recent discussion). The probability of finding a v e at the detector depends on the probability 
of z/ 2 — > v e oscillations inside the earth P 2e . Since in this case the initial state is a mass 
eigenstate z/ 2 , one has a\ = sin 6*0, a 2 = cos8 . In the short pathlength limit, in the case of 
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matter of constant density, one finds from @ 



P2e - {PveUtt = ^2e - sin 2 9 = sin 2 28 (V5-t 2 ) . (16) 

This expression vanishes in the V — > limit because mass eigenstates do not oscillate in 
vacuum. However, it is nontrivial that it is of the order V5-t 2 rather than (V5-t 2 ) 2 ; because 
of this the day-night effect in solar neutrino experiments can be sizeable and must not be 
neglected even when the neutrino pathlengths inside the earth are relatively short. Since a\ 
and d2 are both real, there is no ~ 8-t contribution to the transition probability. 



Another interesting point to notice is that the r.h.s. of ( Tig) is oc 1/E. For sizeable 
pathlengths, the day-night effect increases with neutrino energy when 5 ^> V and decreases 
when 5 <C V (see, e.g., fig. 4 in flOfl). From ( |i~6|) it follows that for short pathlengths it 
always decreases with E, irrespective of the relative magnitudes of V and S. 

Let us consider a few numerical examples. For Vt ~ 84 ~ 1/3, which corresponds to the 
neutrino pathlength inside the earth t ~ 1270 km and 8 ~ V ~ 5.2 x 10~ 14 eV, and assuming 
sin 2 29q ~ 1, the matter-induced oscillation probability (|T6|) is of the order of 10%. For 7 Be 
solar neutrinos (E = 0.862 MeV) the above value of 5 corresponds to Am 2 ~ 1.8 x 10 -7 
eV 2 , which is in the range of the MSW-LOW solution of the solar neutrino problem (for 
a recent analysis of the solar neutrino data see ||11|| ). For the same distance t ~ 1270 km, 
matter effects on the probability of neutrino flavor oscillations would constitute only about 
1%, i.e. an order of magnitude smaller. For typical parameters of the MSW-LMA solution 
and E ~ 10 MeV (typical energy of the 8 B solar neutrinos) eq. (0) yields the probability 
~ 1% for the pathlengths as short as about 100 km, which has to be compared with the 
value 10~ 4 for the earth matter effect on neutrino flavour oscillations over the distance of 
100 km. 

Similar considerations apply to oscillations of supernova neutrinos inside the earth since 
those neutrinos also arrive at the earth as mass eigenstates. If the next supernova explosion 
occurs at such a time that its neutrinos come to a terrestrial detector passing through a 
short distance inside the earth, the earth's matter effects on their oscillations may still be 
quite strong and should be taken into account 0. 

Mass eigenstate neutrinos are also produced in neutral current reactions, i.e. in decays 
of real or virtual Z° bosons; it appears, however, technically rather difficult to produce 
sizeable beams of neutrinos born in neutral current reactions. 



6 Neutrino oscillations due to FCNC 

Neutrinos can oscillate in matter even if their masses are zero or negligible, provided that 
they have flavour changing neutral current (FCNC) interactions |T|. The evolution of the 

4 The importance of matter effects on oscillations of supernova neutrinos inside the earth was recently 
emphasized in |12[; however, the short pathlength limit was not discussed in these papers. 
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system is described by eq. ([!]) with A and B being now both t dependent: 



A(t) = ^fe'N(t), B(t) = V2G F eN(t), (17) 

Here e and e' are the FCNC parameters. In the short-baseline approximation, the transi- 
tion probability P[y a — > 2/&; t) ~ B 2 t 2 . In contrast to the case of the ordinary neutrino 
oscillations in which the off-diagonal element B of the effective Hamiltonian is independent 
of matter density N(t), in the case of oscillations due to FCNC it is proportional to N(t). 
Therefore the matter effect on the oscillation probability (which in this case coincides with 
the probability itself) is of the order {Vt) 2 rather than (VS-t 2 ) 2 . This explains why in this 
case, unlike in the case of the ordinary neutrino flavour oscillations, matter effects can be 
quite sizeable even for as short baselines as that of the K2K experiment, t = 250 km (this 
fact was previously pointed out in [[13], [HJ). The expression P{y a — > v\>\ t) ~ B 2 t 2 also 
explains why in the short baseline limit matter effects do not depend on e' in the leading 
order [O . 



7 Summary and conclusion 

We discussed matter effects in short-pathlength neutrino oscillations and found out the 
reason why these effects generally decrease with decreasing pathlength more rapidly than 
the oscillation effects themselves. This happens because, in the leading order in perturbation 
theory, the transition amplitudes are determined by the off-diagonal terms of the effective 
Hamiltonian which are matter independent. This, however, is not true if the initial neutrino 
state is not a flavour eigenstate, in which case the diagonal terms of the effective Hamiltonian 
also contribute, or if the oscillations are due to flavour changing neutral currents, when the 
off-diagonal terms depend on matter density. 

Initial states which are not flavour eigenstates can be obtained if neutrinos propagate 
significant distances in vacuum before entering matter (or after exiting it). In this case mat- 
ter effects in short-pathlength neutrino oscillations can be strongly enhanced. We discussed 
implications of this observation for atmospheric neutrinos with nearly horizontal trajectories 
and for short pathlength accelerator experiments with the detector placed at large distances 
from the point where neutrinos exit the earth. 

One also deals with initial neutrino states which are not flavour eigenstates when con- 
sidering the earth matter effects on solar neutrinos in the case of the MSW solutions of 
the solar neutrino problem, or oscillations of the supernova neutrinos inside the earth. In 
these cases the initial states are mass eigenstate neutrinos. Matter effects for short neutrino 
pathlengths inside the earth are of the order V5-t 2 rather than (V5-t 2 ) 2 which would be 
expected if the initial state were flavour eigenstate. 

In conclusion, we have shown that matter effects in oscillations of neutrinos traveling 
short distances in matter can be strongly enhanced when the initial neutrino state entering 



S 



the matter is not a flavour eigenstate, or when neutrinos are detected at significant distances 
from the point where they exit the matter. We believe that this is an interesting observation, 
even though the conditions for such an enhancement are generally difficult to realize. A 
notable exception is provided by solar and supernova neutrinos for which these conditions 
are realized as they arrive at the earth in mass rather than in flavour eigenstates. 

The author is grateful to H. Minakata for asking a question during the NOW2000 
workshop which led to the present study, and to A. Yu. Smirnov for useful discussions. 
This work was supported by Fundagao para a Ciencia e a Tecnologia through the grant 
JNICT-CERN/P/FIS/15184/99 and by the TMR network grant ERBFMRX-CT960090 of 
the European Union. 
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